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Let  ft  be  a  finitely  connected  closed  point  set  in  the  conplex  plane 
with  a  piecewise  smooth  boundary  3ft.  The  approximation  of  functions  analytic 
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approximation  at  preselected  nodes  is  discussed.  Attention  is  focussed  on 
simple  methods  for  selecting  an  appropriate  rational  space  and  obtaining  a 
fairly  well-conditioned  rational  basis.  Applications  include  the 
determination  of  conformal  mappings.  Numerical  examples  illustrate  the 
approximation  method. 
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SIGNIFICANCE  AND  EXPLANATION 


Complex  rational  approximation  by  interpolation  has  a  long  history  in  the 

theory  of  approximation.  The  following  numerical  questions,  however,  do  not 

appear  to  have  received  much  attention: 

1)  Given  a  region  on  which  an  analytic  function  shall  be  approximated  by 
rational  functions,  and  given  a  set  of  interpolation  points  on  the 
boundary  of  this  region,  how  should  one  numerically  determine  a  suitable 
rational  space? 

2)  How  does  a  well-conditioned  basis  of  this  rational  space  look? 

3)  If  one  is  free  to  select  interpolation  points  on  the  boundary,  how  should 
they  be  chosen? 

4)  Can  the  selection  of  the  rational  space  be  sinplified  if  one  allows  least 
squares  approximation  instead  of  interpolation? 

The  present  paper  discusses  these  questions. 


ON  COMPLEX  RATIONAL  APPROXIMATION 
BY  INTERPOLATION  AT  PRESELECTED  NODES 

Lothar  Reichel 

1.  Introduction 

Let  ft  be  a  closed  region  of  finite  connectivity  in  the  complex  plane, 
and  assume  that  the  boundary  3ft  is  piecewise  smooth.  Let  f(z)  be  a 
function  analytic  on  ft  and  assume  that  f(z)  is  explicitly  known  on  3ft  or 
on  a  finite  point  set  on  3ft.  The  purpose  of  the  present  paper  is  to  describe 
a  numerical  method  for  determining  a  rational  approximant  r(z)  to  f(z)  on 
ft.  The  method  consists  of  the  following  steps 

1)  select  finitely  many  nodes  z^  in  the  point  set  of  3ft  on  which  f(z) 
is  known. 

2)  choose  a  rational  space  from  which  the  approximant  r(z)  is  to  be 
selected.  The  choice  will  depend  on  the  distribution  of  nodes  z^. 

3)  select  a  well-conditioned  basis  of  the  rational  space. 

4)  compute  r(z)  by  interpolation  or  least  squares  approximation  at  the 
nodes  z^ • 

Complex  approximation  by  interpolation  has  a  long  history  in  the  theory 
of  approximation.  Our  scheme  differs  from  previously  described  methods,  see 
[1],  [3],  [7],  [9],  in  that  we  use  the  selection  of  nodes  as  starting  point. 
This  allows  us  to  treat  cases  when  f(z)  is  known  on  a  finite  point  set  only, 
and  it  also  allows  us  to  let  the  allocation  of  nodes  depend  on  properties  of 
f (z).  In  turn  we  discuss  necessary  and  sufficient  conditions  on  the  selection 
of  rational  space  (section  2),  a  sinple  method  for  choosing  rational  space 
(section  3),  the  selection  of  nodes  (section  4),  choice  of  rational  basis 
(section  5),  and  approximation  on  multiply  connected  regions  (section  6).  We 
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2.  Convergence  results 

Throughout  this  paper  {z  }™  denotes  a  set  of  interpolation  or  least 

Kffll  K*  I 

squares  nodes  on  3Q.  The  set  {w.  .  denotes  a  set  of  not  necessarily 

k  t  n  k  i 

distinct  poles  in  8  ,  the  complement  of  £2  with  respect  to  the  extended 

c 

complex  plane,  and  defines  the  rational  space 

-1  -1  -1  n_1  -1 
(2.1)  Q  :*  span{l,(z-w,  )  ,  <z-w.  )  (z-w  )  *  (z-w  )  } 

n  i#n  i»n  z,  n  x,n 

The  approximation  error  we  measure  in  the  maximum  norm 

I fl,  :*  max  |f (z) | 

3n  zeaa 

The  next  definitions  follow  [9] . 

Definitions. 

Let  the  real  valued  function  c  be  positive  a.e.  on  a  Jordan  arc  y, 

and  assume  /  a(z)|dz|  =  1.  The  direction  of  integration  defines  an 

* 

orientation  on  y,  and  we  let  z  be  the  first  point  of  y.  The  mapping 

F  :  y  +  [0,1],  F(c)  S“  o(c)|dcl  defined  by  integration  along  in  the 

^  n 

positive  direction  has  an  inverse  a.e.  .  For  a  sequence  of  sets  {CfcjJfc. i » 

n  »  1,2,3...  of  points  on  y,  let  for  constants  0  <  d1  <  d2  <  1, 

Nn(d^,d2)  denote  the  number  of  points  of  on  t*ie  8ubarc  of  Y  with 

end  points  F  ^Cdj)  and  F  ^(d2).  If  lim  —  N^td^.d^)  “  d2  -  dj, 

n+» 

Vd^ ,d2,  0  <  d1  <  d2  <  1,  then  the  point  sets  are  said  to  be  uniformly 

distributed  on  y  with  respect  to  o  as  n  ♦  *».  A  set  j=1  P°^n^s 

C jn  on  y  is  said  to  be  eguidistributed  with  respect  to  Y  if 

F  *(£  )  -  p  )  »  — ,  k  ■  2(1 )n.  Sequences  of  equidistributed  sets 

k,n  k-1 ,n  n 

are  uniformly  distributed.  The  definitions  carry  over  to  Jordan  curves  if  we 

•  * 

let  G  be  a  point  on  the  curve,  and  identify  the  curve  with  an  arc  with  t. 


both  as  first  and  last  point 


The  following  theorem  covers  approximation  on  bounded  simply  connected 


regions.  Extensions  are  provided  in  the  remark  below. 
Theorem  2.1 


Let  T  .  and  T  be  piecewise  smooth  Jordan  curves,  T 

nodes  poles  poles 


containing  T  in  its  interior,  and  r  n  T  *  (J.  Let  S  denote 

nodes  poles  nodes 


the  open  region  between  T  ,  and  T  ,  .  Let  0(z)  solve  the-  Mrichlet 

nodes  poles 


problem 


r 


(2.2) 


O(z) 

is 

harmonic  in 

x#: 

y,  z  -  x  +iy. 

0(z) 

is 

continuous  on 

O(z)  =  1 

on 

^  nodes 

0(z)  =  0 

on 

r  , 
poles 

'  u  r 

poles  nodes 


Let  r—  denote  the  outward  normal  derivative  with  respect  to  S.  Then 
d  n 

r  30(  z)  .  .  ,  30(z)  , 

(2.3)  c  :m  I  T" —  Idz  «  -/  —r~ —  |  dz  >  0 


3n 


nodes 


poles 


Let  f(z)  be  a  function  analytic  on  and  interior  to  :=  {z,  0(z)  = 

where  u  is  a  constant  0  <  u  <  1.  Let  {z.  }”.,  n*  1,2,3,...  be  a 

>  1  K ,  n  k=  i 

uniformly  distributed  with  respect  to 


sequence  of  points  sets  on  r 


nodes 


-1  30 

c  r—  as  n  + 


3n 

sets  on  r 


>,  and  let  {w^  n  *  2,3,4,...  be  a  sequence  of  point 


,  n-1 


-1  30 

,  uniformly  distributed  with  respect  to  -c  — .  Let  the  sets 
poles  an 


{w,  )/  ^  define  a  sequence  of  rational  spaces  Q„,  n  *  2,3,4.. 

K,n  k3 I  n 


see 


(2.1).  Then  rn  S  Qn,  uniquely  determined  by  interpolating  f(z)  at  points 
n 


in  the  set  {z.  ,  converges  to  f(z)  on  and  interior  to  T 

k  $  n  K3 1 


for  all 


0  <  P2  <  P,,  as  n  ♦ 


The  rate  of  convergence  is  given  by 


(2.3) 


lim  max 
n+«  zer 


-2n(y  -g  ) 

I f ( z )  -  r  (z)|Vn  <  e  1 

n 


,  n-1 


If  the  sequences  {  z.  }  ,  { w,  }  are  uniformly  distributed  with  respect 

K  |D  )c® *  k  f d  k®  i 

to  another  density  function  such  that 

Vz)  i  j-  -  i  J,  lnls-Wy  _i 1 

k«  1  k*  1 


(2.4a) 

(2.4b) 


lim  V  (z) 
n 


V(z)  is  nonconstant  on  T  , 

nodes 


then  there  is  a  function  g(z)  analytic  on  and  interior  to  r  .  ,  such 

nodes 


that  r  (z)  e  Q  ,  r  (z,  )  *  f(z.  ),  k  ■  1(1)n,  and  Ig-r  I 

n  n  n  k,n  k,n  n  T 


0  as 


nodes 


Proof.  Equations  (2.2)  have  a  unique  solution  U(z)  with  >  0  a.e.  on 
a  9n 

3u 

F node s  and  3n  <  ®  on  I*  Green's  formula  yields  (2.3).  By 

studying  potentials  (2.4a)  Walsh  has  established  the  connection  between  the 

A 

level  curves  of  u  and  the  rate  of  convergence,  see  the  proof  of  Theorem  9  in 

Walsh  [10],  ch.  8.  If  the  limit  potential  V(z)  in  (2.4)  is  nonconstant  on 

T  .  ,  then  there  are  points  C,  e  T  ,  ,  C  6  exterior  of  T  such 

nodes  1  nodes  2  nodes 

that  V(e  )  >  V(c  ).  Let  g(z)  (z-C,)  \  and  let  r  B  Q  interpolate 

g(z)  at  z  -  zkn,  k  «  1 ( 1 )n.  Then  by  [10],  ch.  8, 

fr  .  fJ  .  n  V*nk  "-1  S^nk 

1  k-1  *2  nk  k-1  C1  nk 


and 


Since 


ln|g(S  )  -  r  (C, >  1  -  n(V  (?,)  -  V  (?„))  . 

•  n  i  n  i  n  2 

v  (c  )  -  V  u  )  V(C  )  -  V(c  )  «s  5  >  0,  n  ♦  -  , 

n  i  n  i  l  2 


IgtC,)  -  rn(C , >  I  -  e 


n6 


»,  n 


which  shows  divergence  and  completes  the  proof 


u  r  ) ,  and  c 
nodes  poles  $ 


Remark  2.1 

The  distribution  of  nodes  and  poles  described  in  the  theorem  is  invariant 

under  conformal  mapping,  see  [10],  section  9.12:  Let  $  map  S  conformally 

and  1  to  1  onto  $(S)  and  be  continuous  and  1  to  1  on 

S  T  r  .  .  If  the  node  sets  (z,  },n  „  are  uniformly  distributed  on 

nodes  poles  kn  k=1 

-1  3U 

T  ,  with  respect  to  c  r —  as  n  ",  where  U  solves  (2.2),  then  the 

nodes  3  n 

node  sets  {$(z,  )}”  4  are  uniformly  distributed  on  <MT  .  )  with  respect 

-  kn  k*l  nodes 

“1  ” 

to  c.  t— *-  as  n  ♦  »,  where  U,  solves  the  Dirichlet  problem  analogous  to 
9  on  9 

(2.2)  on  the  mapped  region  $ ( S  u  T 

3U  n-1 

_  ,  r— ®-(z)  |dz| .  Similarly,  if  {w.  },  4  are  uniformly  distributed  on 

?(r.)3n  kn  k**  i 

nodes 

T  ,  with  respect  to  -c  1  as  n  ♦  ®,  then  the  sets  {<9(w,  )},n  ]  are 
poles  3n  kn  k=1 

3U 

uniformly  distributed  on  $(r  0^es>  with  respect  to  -c^  as  n 

Especially  theorem  2.1  holds  also  if  F  ,  is  exterior  to  T  ,  or  if 

nodes  poles 

T  is  a  piecewise  smooth  Jordan  arc. 

nodes 


Remark  2.2 

The  configuration  of  curves  in  theorem  2.1,  may  consist  of  several 

mutually  exterior  curve  pairs  {T  «  ,  T  ,  },  and  the  allocation  of  nodes 

nodes  poles 

and  poles  on  each  pair  can  be  made  independent  of  the  allocation  on  the  other 
pairs.  This  follows  from  the  fact  that  for  each  curve  pair  a  solution  of 
(2.2)  if  extended  to  the  exterior  of  the  curve  pair  would  be  constant  there. 
Its  normal  derivative  on  any  other  curve  would  vanish.  This  remark  follows 
again  from  [10],  ch.  8,  theorem  9  and  its  proof,  which  covers  a  more  general 
situation  than  theorem  2.1  above. 


U 


M 


-6- 


We  close  this  section  by  Indicating  how  theorem  2.1  can  be  used  for 
computing  analytic  continuations.  Let  f(z)  be  known  on  a  curve  ^notjea *  1 

many  physical  problems  one  may  know  that  f(z)  Is  analytic  in  a  specific 
simply  connected  region  B  containing  r  in  its  interior,  lhen  let 

Tpoles  he  boundary  curve  of  B.  The  rational  approximants  computed 

as  described  in  the  theorem  converge  to  f(z)  in  the  interior  of  r  as 

n  ♦  •.  The  next  section  discusses  how  the  nodes  zkn  and  poles  w^n  can  be 
allocated  without  explicitly  solving  the  Dirichlet  problem  (2.2). 


Selection  of  rational  space 


We  discuss  approximation  of  analytic  functions  f(z)  on  simply  connected 

regions.  Our  starting  point  is  the  assumption  that  a  density  function  o  for 

the  interpolation  nodes  { z,  }T  .  on  T  is  known,  and  that  the  nodes 

K ,n  k=l  nodes 

are  eguidistributed  w.r.t.  a.  We  assume  o  >  0  a.e.  on  r  ,  .  If  a  set 

nodes 

of  nodes  {z,  }”  .  is  given  on  T  .  ,  then  we  construct  a  piecewise 

k,n  k=1  nodes 

linear  density  function  such  that  the  nodes  are  equidistributed  w.r.t.  the 
constructed  density  function,  which  we  also  denote  by  o .  A  set  of  poles 
{w  1  defining  the  space  Q  are  obtained  by  solving  (2.2)  as  an  initial 

K  0  n  n —  1  ** 

3U 

value  problem:  U  and  r —  =  a  are  known  on  T  ,  and  we  want  to 

9n  nodes 

determine  other  level  curves  of  U,  on  one  of  which  we  ">llocate  the  poles 


An  initial  value  problem 

Assume  that  T  is  a  smooth  Jordan  curve.  If  T  has  corners 

nodes  nodes 

we  round  them  for  the  present  computations.  If  T  is  an  arc,  we  replace 

nodes 

the  arc  by  a  smooth  circumscribing  curve,  or  we  could  proceed  as  illustrated 

in  section  8.  First  determine  a  set  of  n  -  1  points  { t  }n  1  equi- 

k,n  k“ 1 

distributed  w.r.t.  o.  Let  W  denote  the  conjugate  harmonic  function  to  U 


such  that  w(G  )  =  0.  The  conformal  mapping 
•  ,n 

A 

z  ♦  z  =  $(z)  :=  exp(U(z)  +  iw(z)) 


maps  S,  see  theorem  2.1  for  a  definition,  conformally  on  an  annulus  with 

k-1 

♦  (5.  _)  =  exp(1  +  2iri  - ),  k  =  1(1)n-1.  Now  assume  that  S  is  exterior  to 

^nodes'  conforma^  invariance  noted  in  remark  2.1,  the  poles 

(w  should  be  allocated  so  that  the  points  ^(w,  ),  k  =  1(1)n-1,  are 

K,n  k=i  k,n 

equidistant  on  a  circle  concentric  with  the  unit  circle.  Let  <(>  1  be  the 

*  t+is 

Inverse  of  $ ,  and  let  z  *  e  ,  s,t  e  R.  Then 


(3.1) 


1+2n  i  — 

S.n  =  ♦"  (•  n  )  *  k 


1 ( 1 ) n- 1 


For  a  fixed  t  =  tg,  the  curve  z  =*  z(tg,s),  0  <  s  <  2m,  is  a  level  curve 
of  U.  We  determine  such  level  curves  by  solving  an  initial  value  problem  for 

the  Cauchy-Riemann  equations  for  z  =  z(t,s), 

,,  ,  ,  3z  3z 

3 1  3  s 

Initial  value  problems  for  (3.3a)  are  ill-posed,  but  a  low  accuracy  solution 
suffices  for  our  purpose,  and  generally  we  integrate  few  steps  only.  The  ill- 
posedness  has  not  caused  any  difficulty  in  the  present  application.  On 

F nodes  (3,1)  yields'  with  \  a  *  m  H1)n, 


(3.4) 


3  z  ,  4  . 

17  "'V 


^k+t.n  ^k-1.n 
2As 


Substituting  (3.4)  into  (3.3a),  integrating  in  the  positive  t-direction  with 
Eulers  method  with  At  »  As,  and  denoting  the  computed  approximation  of 


z(1  +  At/S^)  by  yields  the  scheme 


(3.5a) 


Wk,n  s“  Ck,n  *  2  (S+1,n  “  «k-1,n)'  k  "  1(1)n'1'  <n  <0 


«“  C 


n-1 


The  wv  n  He  on  an  approximate  level  curve  of  U,  and  w,  lies 
*  K^n 

approximately  on  the  same  stream  line  as  C  .  Hence, the  w. 

k,n  k,n  are 

3  U 

approximately  equidistributed  w.r.t.  r—  on  a  level  curve  of  U. 

o  n 


Ex.  3.1.  Let 


T  j«  (z(t)  :=  x(t)  +  iy(t),  x(t)  :*  1.75*cos(t) 

node  8 

(3*6)  +  2.625*  cos(2t)  -  2.625,  y(t)  :»  3.0625* sin(t-0. 2) 

+  1.225*  sin(2t)  -  0 .6125*  sin(4t)  +  0.875,  0  <  t  <  2rr }  . 

Allocate  32  points  C.  ,,  equidistantly  w.r.t.  arc  length  on  T  ,  and 

nodes 

integrate  according  to  (3.5a).  Figure  3.1  is  obtained. 


nodes 


Figure  3.1 

The  curve  (3.6)  has  been  used  by  Meiss-Markowitz  [5]  in  a  quite  different 
example. 

■ 

The  integration  (3.5a)  can  be  repeated  by  first  letting  £  :=  w  , 

k,n  k,n 

k  =  1(1)n-1,  and  then  performing  (3.5a).  The  integration  should  be  repeated 

until  the  computed  level  curve  intersects  itself,  and  the  n  should  be 

allocated  on  the  last  non-intersecting  or  near-non-intersecting  computed  level 

curve.  We  motivate  this  by  considering  the  case  where  r  is  the  unit 

nodes 

circle.  Analogous  results  can  be  established  for  more  general  curves 
^ nodes  * 

Ex.  3.2.  Let  T  =  {z  :  |z|  *  1}  ,  and  let  the  interpolation  nodes  be 

nodes 

equidistant  on  T  .  S  is  exterior  to  T  ,  ,  and  level  curves  of  U 

nodes  nodes 


are  circles  I z  t  *  r  >  1.  The  poles  w^  n  will  lie  equidistantly  on  a  circle 
|z|  -  rQ  >  1,  and 

U(z)  -  1  -  ln|z|  •  (ln|zQ|)  1 

«w  M  A* 

At  any  fixed  point  z,  1  <  |z|  <  rQ,  we  have  that  0(z)  increases  with 
|zq|>  Hence,  approximation  of  analytic  functions  on  the  unit  disk  by 
interpolation  at  the  roots  of  unity  is  by  theorem  2.1,  best  done  by  the  family 
of  rational  functions  which  correspond  to  -  00 ,  i.e.  polynomials. 


When  approximation  of  functions  on  the  region  exterior  to  T  .  ,  is 

nodes 


considered,  then  S  is  in  the  interior  of  ^no<jeg*  and  (3.3a)  is  replaced  by 


(3.3b) 


3z  3z 


3t  “  1  3s  * 


The  corresponding  difference  equation  is 


(3-5b>  *  2  'Vi,.  '  k  ‘  '*  'l'  ?0  -  Vi 


Method  (3.5)  as  well  as  other  integration  methods  for  (3.3)  have 
been  studied  in  [6]  for  the  case  when  the  function  to  be  continued 
analytically  is  analytic  in  a  simply  connected  region.  The  analysis  carries 
without  difficulty  over  to  the  present  situation. 

We  conclude  this  section  with  some  computed  examples.  All  computations 
in  this  paper  were  carried  out  on  a  VAX/780  in  double  precision  arithmetic, 
i.e.  with  12  significant  digits. 


Ex.  3.3.  Approximate  f(z)  ;=■  / z-a  on  and  interior  to  the  curve  r  of 

nodes 


example  3.1.,  with  a  : =  -2.6325  +  i  1.425,  see  figure  3.2 


Figure  3.2 


The  branch  of  the  square  root  is  chosen  to  make  f(z)  single  valued  and 

analytic  in  the  plane  cut  along  z  :=  a  +  t,  t  <  0.  For  n  *  1  +  321, 

l  -  1(1)4,  we  allocate  n  interplation  nodes  zy  n  equidistantly  w.r.t.  arc 
length  on  T  ^  .  q32  is  defined  by  the  32  poles  wjf33  on  fi9ure  3.1. 

°33+32*'  1  =  1,2,3  arC  deflned  by  letting  Wj/33+32jt  !~  w(j  mod  32)+1,33  Vj* 

Let  rn(z)  denote  the  element  in  Qn  such  that  rn(zk#n)  *  ^zk,n>' 

k  =  1 ( 1 )n.  In  figure  3.3  the  computed  errors  are  marked  with  dots. 

We  note,  in  passing,  that  the  nodes  of  course  do  not  have  to  be  allocated 
exactly  equidistantly  w.r.t.  arc  lengths.  Nodes  could  be  marked  sufficiently 
accurately  with  a  light  pen,  and  also  equations  (3.3a,b)  are  sufficiently 
simple  to  allow  an  approximate  graphic  solution. 


I 
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Figure  3.3 


Ex.  3.4.  In  theory,  it  is  also  possible  to  approximate  the  function  f(z)  of 

A 

example  3.3  by  polynomials.  Let  {*.  }  denote  a  set  of  n  Fejer  points 

k,n  k«i 

for  T  _  .  For  their  definition,  see  [3J  or  [9] .  Figure  3.4  shows  the 

nodes 

points  (z!  . marked  with  crosses  on  T  _  Interpolation  of  f(z) 
k,  120  k-1  nodes 

in  n  points  defines  a  polynomial  pn(z)  of  degree  <  n,  and  the  polynomial 

sequence  {p  (z)}  converges  maximally  to  f(z),  i.e.  p  (z)  converges 
n  i  n 

exponentially  to  f(z),  n  ♦  ",  and  there  is  no  polynomial  sequence  with  a 
higher  exponential  rate  of  convergence  to  f(z). 


3- 


Figure  3.4 


Computation  of  pn  for  some  n  gave  the  following  table 

n  lf-p  I 

nodes 

40  0.60 

80  0.42 

120  0.42 

160  impossible  to  evaluate. 

A  slow  rate  of  convergence  is  combined  with  difficulties  of  accurately 

evaluating  pn(z)  for  large  n.  A  Lagrange  polynomial  basis  was  used,  and 
*  F 

the  Fejer  points  zfc,n  were  determined  with  4  significant  digits. 

I 

Ex.  3.5.  Approximate  f(z)  :*  -  / <z-z, ) (z-z,)  on  and  exterior  to  T 

z  i  2  nodes 

defined  by  (3.6),  where  z.  :=  -  4-2i,  z,  =  i,  see  figure  3.5. 


Figure  3.5 


Figure  3.6 


The  branch  of  the  square  root  is  chosen  so  that  zf(z)  is  analytic  in 

the  finite  plane  cut  between  z1  and  z2*  We  wish  to  approximate  f(z)  by 

interpolation  in  33  nodes  equidistant  w.r.t.  arc  length.  Figure  3.6  shows  32 

points  C  allocated  equidistantly  w.r.t.  arc  length  on  T  ,  and  also 

j  t 33  nod€8 

32  poles  wjt33  obtained  by  applying  (3.5b)  once,  r33  denotes  the  rational 

interpolant  to  f(z)  in  Q33,  see  below. 

Allocate  65  Interpolation  nodes  on  T  equidistantly  w.r.t.  arc 

nodes 

length.  Let  wj+32  65  !“  wj  33'  3  "  1(1)32.  This  defines 

Let  r^5  be  the  interpolant  to  f(z)  in  Q^.  We  obtain 

n  I  f-r  1  _ 

n  nodes 


-QLZ- 

4.  Selection  of  nodes 

If  f(z),  the  function  to  approximate,  is  known  on  a  discrete  point  set 

only,  there  may  be  no  choice  to  make.  In  this  section  we  assume  that  f(z) 

is  known  everywhere  on  T  .  .  If  no  knowledge  of  the  location  of  the 

nodes 

singularities  of  f(z)  is  available,  we  want  an  allocation  of  nodes  such  that 

the  orthogonal  distance  from  a  point  of  T  to  level  curves 

nodes 

{z  :  u(z))  =  p} ,  p  >  0,  is  approximately  constant  for  all  points  of  F 

nodes 

Then  the  rate  of  convergence  will  depend  on  the  distance  from  T  to  a 

nodes 

singularity  of  f(z)  closest  to  T  This  is,  to  a  first  approximation, 

nodes 

achieved  by  allocating  the  nodes  equidistantly  w.r.t.  arc  length. 

Knowledge  about  the  location  of  the  singularities  of  f(z)  can  be  used 

by  allocating  more  nodes  on  parts  of  T  close  to  a  singularity.  The 

nodes 

level  curves  of  U(z)  will  be  close  to  T  where  the  node  density  is 

nodes 


highest.  Example  7.2  provides  an  illustration. 


5.  Rational  basis 


The  basis  implicit  in  the  definition  (2.1)  of  Qn  is  generally  ill* 
conditioned.  If  the  nodes  and  P«l®*  wk,n  are  naar-equidistributed 

with  respect  to  |c  on  rnode8  rpoles'  respectively,  then  the 


n-1  z-z 


k.n 


basis  t0(s)  -  1,  t^z)  s-  n  z  _z 


n-1  z,  -w, 
n  _k*I! — Lifl 
z-w. 


,  j  -  1 ( 1 ) n- 1 ,  is 


k«1  j,n  k,n  k*1  k,n 
k<j 


fairly  well-conditioned.  A  condition  number  of  a  basis  we  define  following 

n-1 

Gautschi  [  ],  with  the  map  F  s  h”  ♦  Q  t  a  ♦  £  a  £  (*),  where  a  * 

n  n  “  ic  K  * 


k-0 


n  — w 

( aQ ,a1 , . . .  ,«n_  1 ) .  Equip  vf1  with  the  maximum  norm  I  j»l.  :*  max|a  |  and 

.  The  induced  operator  norms  are 


Qn  with  the  norm  I  Ij. 


nodes 


n-1  n-1 

IF  I  max  I  £  a  £(*)•_  ■  max  ][  |£  (z)|  , 

n  I  a*  *1  0  nodes  z€T  .  k«0 

—  •  nodes 


IF  ^1  i-  max 

n  _1 

'1  W*"r  „  •< 

0  nodes 


111.  <  2  , 


and  the  condition  number  of  the  basis  £^(z)  is 


cond  F  i«  IF  I  IF  I  <  2  max  £  |£  (z)|  <  2n  •  max  It  (z)l_ 

^  n  n  zer  k”0  0<k<n-1  nodes 

nodes 


n-1 


Assuming  that  the  z^,  w^  are  uniformly  distributed  with  respect  to 
-1  30. 

£  have  for  k  0, 

1  n_1 

ln|£  (z)|  -  (n-1)(—  (  l  lnlz-z  |  -  ln|z  -z  |) 
k  n-1  jm1  j,n  k,n  j,n 


1  n_1 

-  -  ^Vn-^.n-l1 


)) 


<n-1)  (/rnodeslnlz-«lc’1|f”  <«  >  I  ld<i  “  Jrnodeslnl*H,n-< I0"'  <Ol  |dd  + 


-1 ,3U 


-17- 


-  /  j,  In  |  z-C  |c 

poles 


1511 


'*l+/r  ,  lnlsk,n-t|c 

poles 


3n 


for  any  fixed  z  e  Y  ,  z  ?  z  as  n  +  •.  Hence,  | 

nodes 

0(1)n-1,  for  any  fixed  z  €  Y  ,n+®.  This  shows  that 

nodes 

{A.}"  1  is  reasonably  well-conditioned. 


6.  Multiply  connected  regions 

We  consider  the  approximation  of  a  function  f(z)  on  an  exterior  doubly 
connected  region.  Generalization  to  more  complicated  regions  are  immediate. 
Let  f(z)  be  analytic  on  and  exterior  to  the  piecewise  smooth  curves  y^  and 
Y2.  Y ±  bounds  the  region  Q^,  see  figure  6.1. 


Figure  6.1 


We  wish  to  approximate  f(z)  by  a  rational  function  which  interpolates 

f(z)  in  nodes  on  Y1  and  Y2«  During  the  allocation  of  poles  we  can  regard 

the  problem  as  if  it  were  composed  of  the  two  simpler  subproblems: 

Approximate  f^fz),  analytic  on  and  exterior  to  Y^»  1  •  1,2,  and 

f  («)  «  o.  Consider  the  case  i  *  1.  Assume  n  nodes  {z,^}T  .  and  a 
2  k  k*1 

density  function  are  known  on  Y 1  •  The  method  of  section  3,  yields 

poles  w^1^,  k  =  1(1)n-1,  in  and  this  defines  the  space 

(6.1)  Q(1)  :»  span{l,  (  z-w! 1 }  j"1 , . . . ,  II  (z-w*1’)'1} 

n  1  .  ,  k 

k“1 

and  basis 


(6.2) 


r 


*:”<■> 


*"’<*> ■■ 


n-1 

IT 


z-z 


(1) 


(1)  (1) 


j-1  zu  -w 

jA 


,  (1)  (D 

n-1  z.  -w, 

JS _ ] _ 

(1) 


j 


n 

3-i 


z— w , 


1  ( 1 ) n- 1 
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(2) 


Repeat  for  problem  i  =  2.  Nodes  ,  k  =  1...m  are  assumed  to  be  known  on 


Since  =  0,  the  rational  space  we  are  to  construct  does  not  have 


to  contain  constants-  We  can  identify  the  nodes  z 


(2) 


with  the  points  C, 


(2) 


of  section  3  and  obtain  poles  ,  k  =  1,...,m,  defining  the  space 


m 


(6.3) 
and  basis 


Q(2)  :=  spanUz-wJ25)-1,...,  II  <z-w<2>)~1} 
m  1  k=1  k 


ll2){z)  :1 

k 


m 

n 


(2) 


_L 


.  ,  (2)  (2)  .  . 

£  Zk  “Zk  3=1 
jA 


(2)  (2) 
m  z  -w . 

n  — - rih—  ,  j  =  i  ( i  )m 


(2) 

z-w  . 

D 


To  solve  the  original  problems,  we  select  the  function  r  e  Q  ^ 1  *  ®  c/21 

n,m  n  m 

which  interpolates  f(z)  in  u  Convergence  results 


analogous  to  theorem  2.1  can  be  shown,  see  remark  2.2,  provided  that  both 


n,m  +  “.  The  basis  (Z 


(1) 


0 


.(1)  .  ( 2)  (2) .  .  ,  .  . 

,Z  ,,  l,  }  is  fairly  well- 

n—  I  i  m 


( 1 ) , 


conditioned  on  Y  lJ  y  ,  under  the  assumption  that  H l  I  ,  k  =1(1)n-1, 

'  *  k  Y, 

(2) 

and  IJ l  II  ,  k  =  1(1)m  are  small.  The  method  of  proof  is  similar  to  that 

k  Yt 


used  in  section  5.  The  assumption  is  reasonable  due  to  the  relation  between 
nodes  and  poles. 


Ex.  6.1.  Let  f(z 


)  =  /{ lz+  ~)(z+  7)  •  (z+2)  +  / ( z-2-  7) 


2, (z-2+  j)  •  (z-2)"1 


2  2 

where  the  branches  are  selected  so  that  the  first  term  is  analytic  in  the 


5  3 

complex  plane  cut  along  the  line  segment  between  z  =  -  —  and  z  =  -  — ,  and 


the  second  term  is  analytic  in  the  plane  cut  along  the  line  segment  between 


z  =  2  ♦  and  z  =  2  -  Approximate  f(z)  in  the  region  exterior  to  both 


curves 


Y1  =  (z  =  2  +  cos(t)  +  i*sin(t)(2  +  cos(t))  ,  0  <  t  <  2n} 


y^  =  [z  =  -2  -  sin(t)(2  +  cos(t))  +  i»cos(t),  0  <  t  <  2v } 


Figure  6.1  shows  Y  ^ ,  Y2  and  the  branch  points  of  f(z)  marked  with  crosses 


inside  Y 1  and  y 


M 


M 
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mm 


Figure  6.2 


He  wish  to  interpolate  f(z)  in  31  nodes  z^  equidistant  w.r.t.  arc 

(2) 

length  on  and  in  30  nodes  zfc  equidistant  w.r.t.  arc  length  on  Y2* 

Figure  6.2  shows  30  points  equidistant  w.r.t.  arc  length  on  Y  ^  #  the 

(2)  (2)  { 1 )  (2) 

30  points  zfc  “  on  Y2»  and  the  poles  w^  and  w^  obtained  by 

applying  (3.5b)  twice  for  each  curve  y^.  Denote  the  interpolating  rational 

approxieant  by  r(z).  Then 


The  close  connection  between  rational  approximation  and  approximation  by 
rational  harmonics,  see  [11],  suggests  applications  to  the  numerical  solution 
of  Dirichlet  problems  for  the  Laplace  equation  on  simply  and  multiply 
connected  regions.  We  discuss  in  some  detail  the  special  Dirichlet  problems 
whose  solution  yields  a  conformal  mapping  from  a  simply  connected  region  to 
|w|  >  1  or  to  | w 1  <  1. 

Let  ft  be  simply  connected  with  boundary  3ft  =  T  and  complement 

nodes 

ft  .  Allocate  nodes  z.,  j  =  1(1)2n-1,  on  T  for  some  n.  Determine 

c  3  nodes 

n-1  points  C  on  T  ,  such  that  the  distribution  functions  for  z. 

j  nodes  3 

and  c;  .  agree.  Allocate  n-1  poles  w.  in  ft  by  application  of  (3.5a)  or 
3  3  c 

(3.5b).  This  defines  the  2n-1  -  dimensional  space  of  harmonic  functions 

(7.1)  span{ l  (z) ,  Re(£-(z)),  Im(fc  (z) ) , . . . ,Re( i  (z)),Im(Z  ,(z))}  , 

u  1  I  n- 1  n- 1 

where 


(7.2) 


V*  : 


=  1 


l . (z)  :  = 

3 


n-1  z-z  n-1  z  -w, 

_  2k  _  2k  k  . 

ft  -  IT  -  ,  -j  =  1(1)  n—  1 

k=1  zj  2k  k= 1  t 


An  approximate  solution  to 


is  obtained  by  solving 


j  Au  =  0  in  ft 
\  u  =  g  on  3ft 


n-1 

(7.3)  Re(  l  a.ll,  (z.))  =  g(z.)  ,  j  =  1(1)2n-1 

k=0  k  k  3  3 

for  the  a^.  This  system  may  not  always  be  uniquely  solvable,  which  is  why. 


in  the  computed  examples  below,  we  solved  (7.3)  by  singular  value 


decomposition  of  the  matrix.  In  none  of  the  computed  examples  the  matrix  was 
near-singular . 

The  special  choice 

(7.4)  g(z)  =  -ln|z| 

leads  to  approximate  conformal  mappings  for  £2  provided  that  0  ^  3  £2. 


Ex.  7.1.  Let  £1  be  the  bounded  region  of  figure  3.5.  Compute  a  conformal 

mapping  <f>  :  £2  >  |w|  <  1,  such  that  >(0)  *  0.  Allocate  129  nodes  z^  on 

32 

3£2  equidistantly  w.r.t.  arc  length.  Use  the  poles  shown  in  figure 

3.1,  and  define  wj+32  !“  wj*  3  “  1(1)32.  Thus  we  obtain  a  basis  (7.2)  for 
n  =  65,  and  solve  (7.3)  with  g(z)  defined  by  (7.4).  11118  yields 

64 

♦,_(x)  :*  z  exp(  l  a  £  (z)) 

03  ,  -  K  k 

k=0 

which  approximates  $.  Figure  7.1  shows  $ecO£2)  and  a  reference  circle  of 
radius  1 . 
resolution 


63 


1.  The  error  l|$65(z)|  -  1lgn 


7*10 


-5 


is  well  below  the 


of  the  picture. 


Figure  7.1 


Ex.  7.2.  Let  f!  be  the  exterior  of  the  curve  in  figure  7.2.  Compute  a 


conformal  mapping  tj/  :  ft  +  |w|  >  1,  such  that  ^(°°)  =  °°.  Allocate  65 
Zj  on  3ft,  equidistantly  w.r.t.  arc  length,  and  use  the  32  poles  w^ 
figure  3.6.  Solving  (7.3),  with  g  defined  by  (7.4),  one  obtains  the 
approximate  map 


n- 1 

(7.5)  ^  (z)  z  exp(  £  a  Z  (z)) 

n  .  .  K  k 

k=0 


for  n  =  33. 


nodes 

of 


Figure  7.2  Figure  7.3 

Figure  7.3  shows  ^^(Sft)  and  a  reference  circle  of  radius  1.1. 

-2 

H  I4i  (z)  |  -  IS  =3*10  .  >|>  (3ft)  intersects  itself  at  the  blob,  which  is 

33  Hu  33 

roughly  the  image  of  the  part  of  3ft,  which  lies  strictly  interior  to  the 

convex  hull  of  3ft.  This  suqgests  that  in  order  to  achieve  higher  accuracy 

n  3ti 

more  nodes  should  be  allocated  between  the  points  z f  — ' )  and  z(~~)  in  figure 


7.1.  The  allocation  of  nodes  to  be  described  next  is  simplified  by  the  fact 


that  the  parameter  t  in  (3.6)  satisfies 


d(arc  lenth) 


constant. 


0  <  t  <  2ir.  (Figure  2.2  shows  32  points  on  3ft  equidistant  w.r.t.  t.)  Thi 
simplifies  the  determination  of  the  nodes  and  poles  to  be  used,  but  is  not 
essential  for  the  discussion,  why  we  chose  not  to  use  this  fact  in  the  first 
part  of  this  example.  Figue  7.1  shows  32  points  marked  with  crosses  on  T 
and  allocated  equidistantly  with  respect  to  the  boundary  parameter  t,  see 


(3.6).  In  figure  7.4  we  have  allocated  11  points  equidistantly  with 

ir 

respect  to  t  for  0  <  t  <  ~,  44  points  equi -distantly  with  respect 

ir  3* 

to  t  for  —  <  t  <  and  ^ 1  points  equidistantly  with  respect  to 

t,  “  <  t  <  2ir .  By  (3.5b),  we  obtain  66  poles  Wj ,  and  finally  we  allocate 

133  nodes  on  T  having  the  same  distribution  as  the  This  yields  the 

mapping  ♦,-(*)•  Figure  7.5  shows  $  (D  and  a  circumscribed  concentric 
o7  67 

reference  circle.  I|tg7(z)|  “  Up  =  3*  10  3. 


Figure  7.4 


Figure  7.5 


The  abbreviation  error  is  near  the  resolution  of  the  plotter.  Double  the 


number  of  poles  in  figure  7.4  by  Wj+6g  :=  Wj,  j  =  1(1)66.  This  yields 


4  x  66  +  1  =  265  harmonic  basis  function,  which  we  determine  by  interpolation 


at  265  points  z^,  which  we  allocated  with  the  same  density  function  as  we 


-6 


used  for  (J>67.  This  gives  $133(z)  and  ,l$133(z)l  -  Up  =  6*10  .  Figure 


7.6  shows  ^33^^  and  a  reference  circle  of  radius  1.1. 


Figure  7.6 


8.  Least  squares  approximation 

We  consider  an  example,  where  the  distribution  of  nodes  and  poles  does 
not  agree  very  well  with  the  conditions  of  theorem  2.1,  and  point  out  that  in 
such  cases  least  squares  approximation  may  give  higher  accuracy  than  inter** 
polation. 

r 2 

Approximate  f(z)  :  =*  v  z  - 1  on  3(1  **  [-5i,5i],  by  using  function  values 

k-1 

at  equidistant  nodes  z  s“  i(10  -  -  5),  k  ■  1(1)m,  on  30.  Hie  branch 

K|DI  IB 

of  the  square  root  is  chosen  to  make  f(z)  analytic  in  the  finite  plane  cut 
on  the  real  axis  from  1  to  <*  and  from  -1  to  -«*.  Hie  poles  we  allocate 
in  a  simple  manner:  for  n  even,  let  for  some  s  >  0 


f  W2k,n  !"  1(10  iSf  ‘  5)  +  8  ' 

/  k  -  1(1)  |  . 

\W2k-1,n  ’-A<1o2Sif-5>  ' 

First,  consider  the  selection  of  s.  For  s  large  the  distribution  of 
nodes  z^m  and  poles  w^  „  does  not  agree  at  all  with  the  distribution 
suggested  in  theorem  2.1.  In  fact,  for  s  ■  *,  i.e.  polynomial 
approximation,  approximation  by  interpolation  diverges.  This  follows  from  the 
similarities  of  our  approximation  problem  and  the  classical  example  of  Runge, 
see  [2] .  On  the  other  hand,  for  s  >  0  small  the  rate  of  convergence  becomes 
unnecessarily  slow. 


Ex.  8.1.  Let  n  =  m  +  1  and  compute  the  rational  approximant  rn(z)  to 
f(z)  for  s  •  1,2,3. 
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Figure  8. 1 

Only  points  with  markers  in  figure  6.1  correspond  to  computed  approximation 
errors.  g 

Example  8.1  shows  that  the  selection  of  s  is  of  some  importance,  c.f. 
also  example  3.2.  When  it  is  difficult  to  determine  an  appropriate  allocation 
of  poles,  a  crude  determination  of  poles  combined  with  discrete  least  squares 
approximation  at  the  nodes  Zy  m  can  be  a  good  strategy.  For  a  justification 
when  the  poles  all  are  at  00 ,  see  [8].  Approximation  by  rationals  with  a 
different  but  fixed  distribution  of  poles  can  be  treated  similarly  as  in 
[8] .  Given  m  nodes  Zy  m  and  a  distribution  of  poles  w^  n,  one  generally 
does  not  know  a  priori  how  to  select  the  ratio  m/n.  One  has  to  select 
several  values  of  n  and  select  the  best  of  the  computed  approximants,  see 
[8] .  The  difficulties  in  choosing  m/n  are  illustrated  in  the  next  example. 

Ex.  8.2.  Consider  the  same  approximation  problem  as  in  example  8.1  with 
s  =  2,  but  let  "*  =  0.9.  For  m  =  25,  51,  25,  we  let  n  be  the  even 
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integer  closest  to  0.9m.  The  in  figure  8.2  show  the  approximation 


error.  Ttie  dots  correspond  to  interpolation  n  *  m-1  and  are  the  same  as  in 
figure  8.1. 


Figure  8.2 


Acknowledgement 

I  wish  to  thank  Fred  Sauer  for  help  with  the  handling  of  computers  and 


programs 


REFERENCES 


Curtiss,  J.  M.,  Convergence  of  complex  Lagrange  interpolation  on  the 

locus  of  the  interpolation  points,  Duke  Math.  J.,  32,  187-204,  1965. 

Dahlquist,  G.  and  Bjorck,  A.,  Numerical  methods,  Prentice  Hall,  Englewood 
Cliffs,  NJ,  1974. 

Gaier,  D. ,  Vorlesungen  uber  Approximation  im  Komplexen,  Birkhauser,  Basel 
1980. 

Gautschi,  W. ,  Questions  of  numerical  condition  related  to  polynomials,  in 
Recent  advances  in  numerical  analysis,  eds.  de  Boor,  C.  and 
Golub,  G.  H.,  Academic  Press,  New  York  1978. 

Meiss,  Th.  and  Markowitz,  U.,  Numerische  Behandlung  partieller 
Differentialgleichungen,  Springer,  Berlin,  1978. 

Reichel,  L.,  A  numerical  method  for  analytic  continuation  of  conformal 
mappings.  Report  TRITA-NA-8 118,  Dept,  of  Comp.  Sci.,  Royal  Institute 
of  Technology,  Stockholm,  Sweden,  1981. 

Reichel,  L.,  Some  computational  aspects  on  rational  approximation  in  the 
complex  plane.  Report  TRITA-NA-8 1 1 1 ,  Dept,  of  Comp.  Sci.,  Royal 
Institute  of  Technology,  Stockholm,  Sweden,  1981. 

Reichel,  L.,  On  polynomial  approximation  in  the  uniform  norm  by  the 
discrete  least  squares  method,  MRC  Technical  Summary  Report  #2472, 
Mathematics  Reseach  Center,  University  of  Wisconsin-Madison,  1982. 

Smirnov,  V.  I.,  and  Lebedev,  N.  A.,  Functions  of  a  complex  variable, 
Iliffe  Books,  London,  1968. 

Walsh,  J.  L. ,  Interpolation  and  approximation  by  rational  functions  in 
the  complex  domain,  American  Mathematical  Society,  Providence,  R.I., 
1935. 


-3 


I 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  f»k«n  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


«.  REPORT  NUMBER 

#2514 


4.  TITLE  (and  Submit) 


s\v> 


On  Complex  Rational  Approximation  by 
Interpolation  at  Preselected  Nodes 


7.  AUTHORf*) 


READ  INSTRUCTIONS 
BEPORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


S.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Summary  Report  -  no  specific 
reporting  period 


S.  PERFORMING  ORG.  REPORT  NUMBER 


•.  CONTRACT  OR  GRANT  NUMBERftJ 


Lothar  Reichel 


DAAG29-80-C-0041 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Mathematics  Research  Center,  University  of 
610  Walnut  Street  Wisconsin 

Madison.  Wisconsin  53706 


II.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

U.  S.  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park,  North  Carolina  27709 


4.  MONITORING  AGENCY  NAME  4  ADDRESS (it  dlllerent  from  Controlling  Olliea) 


•  0.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  *  WORK  UNIT  NUMBERS 

Work  Unit  Number  3  - 
Numerical  Analysis  and 
Scientific  Computing 


12.  report  oate 
April  1983 


is.  number  of  pages 

31 


IS.  SECURITY  CLASS,  (ol  thlo  rope rt) 

UNCLASSIFIED 


16.  DISTRIBUTION  STATEMENT  (ol  thlo  Report) 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  the  abo tract  entered  In  Block  20.  It  dl Heron t  horn  Report) 


19.  KEY  WORDS  (Continue  on  reverse  aide  it  necessary  and  Identity  by  block  number) 

rational  approximation,  interpolation,  least  squares  approximation, 
conformal  mapping,  analytic  continuation. 


20.  ABSTRACT  (Continue  on  reverse  side  it  nscessory  end  Identify  by  bloc*  number) 

Let  ft  be  a  finitely  connected  closed  point  set  in  the  complex  plane  with 
a  piecewise  smooth  boundary  3ft.  The  approximation  of  functions  analytic  on  ft 
by  rational  functions  determined  by  interpolation  or  least  squares  approximation 
at  preselected  nodes  is  discussed.  Attention  is  focussed  on  simple  methods  for 
selecting  an  appropriate  rational  space  and  obtaining  a  fairly  well-conditioned 
rational  basis.  Applications  include  the  determination  of  conformal  mappings. 
Numerical  examples  illustrate  the  approximation  method. 


nn  FORM 

W  |  JAN  73 


EDITION  OF  I  NOV  46  IS  OBSOLETE 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


